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Abstract 

We consider evolutionary equations of the form ut — F{u, w) where w = D~ 1 D y u is the 
nonlocality, and the right hand side F is polynomial in the derivatives of u and w. The recent 
paper [5j provides a complete list of integrable third order equations of this kind. Here we 
extend the classification to fifth order equations. Besides the known examples of Kadomtsev- 
Petviashvili (KP), Veselov- Novikov (VN) and Harry Dym (HD) equations, as well as fifth 
order analogues and modifications thereof, our list contains a number of equations which 
are apparently new. We conjecture that our examples exhaust the list of scalar polynomial 
integrable equations with the nonlocality w. The classification procedure consists of two 
steps. First, we classify quasilinear systems which may (potentially) occur as dispersionless 
limits of integrable scalar evolutionary equations. After that we reconstruct dispersive terms 
based on the requirement of the inheritance of hydrodynamic reductions of the dispersionless 
limit by the full dispersive equation. 

MSC: 35L40, 35Q5f , 35Q58, 37K10, 37K55. 

Keywords: dispersionless equations, hydrodynamic reductions, dispersive deformations, 
integrability. 
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1 Introduction 



The classification of integrable 1 + 1 dimensional scalar evolutionary equations, 

u t = F(u), 

has been (and still is) a subject of active research within the soliton community. Here u(x, t) 
is a scalar potential, and F denotes a differential expression which depends on x-derivatives 
of u up to some finite order. Although the general classification problem is still out of reach, 
quite a few important results were obtained under various additional assumptions on F (such 
as polynomiality, linearity in the highest derivative, etc). We refer to the review article [12] 
for a detailed discussion of the classification techniques involved, extensive lists of integrable 
equations within particularly interesting subclasses, and references. 

In this paper we apply the novel perturbative approach outlined in [5j|6] to a similar problem 
in 2 + 1 dimensions, the area where very few classification results are currently available. The 
main challenge of higher dimensions is the non-locality of scalar evolutionary integrable equa- 
tions: the corresponding right hand side F must contain nonlocal variables whose differential 
structure was clarified in [13]. Here we consider equations of the form 

u t = F(u,w) (1) 

where u(x, y, t) is a scalar field and w = D~ l D y u is the simplest nonlocality (equivalently, w can 
be introduced via the relation w x = u y ). We assume that the right hand side F is polynomial 
in the x- and y-derivatives of u and w, while the dependence on u and w themselves is allowed 
to be arbitrary. The paper [6] provides a complete list of integrable third order equations of the 
form ([1]), 

u t = ipu x + ipuy + r]w y + e(...) + e 2 (...), (2) 

where tp, ip and r\ are functions of u and w, while the terms at e and e 2 are assumed to be 
homogeneous differential polynomials of the order two and three in the derivatives of u and 
w (one can show that all terms at e have to vanish). We use the following weighting scheme: 
u and w are assumed to have order zero, their derivatives u x ,u y ,w x ,w y are of order one, the 
expressions u xx , u xy , u yy ,w yy ,u x ,u x u y , u y ,u x w y , u y w y ,w y are of order two, and so on. Assuming 
that the dispersionless limit of the equation ([2]), 

U t = (PU X + lf)Uy + TJWy, W X = Uy, (3) 

is linearly nondegenerate (the property to be clarified in Sect. 2.2), and satisfies the condition 
7] 7^ (which is equivalent to the requirement that the dispersion relation of the system (|3|) 
defines an irreducible conic), we have the following result: 

Theorem 1 JBjj Up to invertible transformations, the examples below provide a complete list of 
integrable third order equations (0) with r] ^ whose dispersionless limit is linearly nondegen- 
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erate: 

KP equation 
modified KP equation 

Gardner equation 
VN equation 

modified VN equation 

HD equation 

deformed HD equation 
Equation E\ 
Equation E^ 



u t = uu x + Wy+e u xxx , 

u t = (w - u 2 /2)u x + w y +e 2 u xxx , 

u 2 + 5u)u x + w y +e 2 u xxx , 



u t = {J3w 



P 2 



UW)y-\-€ Uyyy j 



3< 



U t = (Uw)y+e 2 Uyy ~ - — 



u t 



~2WUy + UWy- 



4 u 



,2 ^ 



u \u 



6 e 2 (\ 

—^u x — 2wu v + uw v — 

u u \u 

32„,2\ 



u t = ((3w + l u 2 )u x - 3(3uu y + w y + e z [B*(u) - (3u x B 2 (u)\, 
4 
3 



f3 2 u 3 u x + (w- 3f3u 2 )u y + uw y + e 2 [B 3 (u) - (3u x B 2 (u)], 



here B = f3uD x — D v , f3=const. 



The main result of this paper is a generalisation of the above classification to fifth order 
equations, 

u t = vu x + ^u y + nw y + e(...) + e 2 (...) + e 3 (...) + e 4 (...), (4) 

where the terms at e k are assumed to be homogeneous differential polynomials of the order k + 1 
in the derivatives of u and w, respectively. We also assume that the e 4 term depends on at least 
one of the possible fifth derivatives u xxxxx , u XXXX y, . . ., i. e. the equation (d|) is of order 5. 

Theorem 2 Up to invertible transformations, the examples below provide a complete list of in- 
tegrable fifth order equations ^ with n ^ whose dispersionless limit is linearly nondegenerate: 



BKP equation 
CKP equation 

HD§ equation 

Equation E% 

Equation E± 



u t = 5(u 2 + w)u x + buuy - bw y + 5e 2 (uu xxx + u xxy + u x u xx ) + e u xxxxx , 

5 

u t = 5(u 2 + w)u x + buuy - 5w y + 5e 2 (uu xxx + w xxx + -u x u xx ) + e 4 u xxxxx , 



Ut 
Ut 

+ 5< J 



lbwuy — 5uw y + 5e 



u xxy 3 / U X Uy 

u 2 



U \ M z 



2 

: '1 



2u 2 



ii- 



Ari 1 — %■ + 5(3w — \)u v — 5uw v 

■7/0 u 



7 



2u 2 \ u 



-I + 



3 / U X Uy 



U \ U" J x 



2u 2 



u- 



ut = 47 2 ^ + 5(3w — \)u y — 5uw y 



+ 5e 2 



u° 

' l. (L 

3u \u 3 

— (- 

2u 2 { u 2 



u x 

7^7 

u' 



+ 



xxy 

15 ( (2uu xx -3u 2 x ) 2 

16 V u 8 



4u 4 



\2UUn 



3u 2 x ) 
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We point out that the last two examples from Theorem 2 are apparently new. The equation 
E% can be viewed as a deformation of the fifth order Harry Dym equation HD5: it reduces 
to HD5 when 7 = 0. Although each equation appearing in Theorems 1-2 gives rise to an 
integrable hierarchy, the corresponding higher flows will not belong to the class ([1]): they will 
necessarily have a more complicated nonlocality. Preliminary calculations suggest that there 
exist no seventh order equations of the form ([1]). This leads to the following 

Conjecture Up to invertible transformations, Theorems 1-2 provide a complete list of integrable 
evolutionary equations of the form (CP which are polynomial in the derivatives of u and w. 

Remark. The assumption of polynomiality is essential: there exist examples of integrable 
equations of the form (pQ) where the right hand side F is an infinite series in e. As an illustration, 
let us consider integrable differential-difference equations of the Toda lattice, 

v t = vA_(w), w x = A + (v), 

where 

A_(w) = w (y) - ™(y - e ) ^ A+( „ ) = "fo + *)-«(*). 

Introducing the variable u by the formula A + (v) = u y , one can rewrite the equations of the 
Toda lattice in such a way that the nonlocality w will be of the required form, 

u t = D y X A + (A+ 1 (u y )A^(w)) , w x = u y . 

Expanding the first equation in powers of e one obtains an infinite series, 

U t = UWy + —(UWyy)y + W X =Uy. 

Examples of this type will be outside the scope of this paper. 

The structure of the paper is as follows. Following [6] , in Sect. 2.1 we review the classification 
of integrable quasilinear systems of the form ([3]). In Sect 2.2 we outline the general procedure 
which, starting with an integrable dispersionless system, allows one to systematically reconstruct 
dispersive corrections. This procedure is applied in Sect. 2.3 to the case of fifth order equations 
dU). For the reader's convenience, in Sect. 3 we present Lax pairs for all equations appearing in 
Theorems 1-2. 

2 Proof of Theorem 2 

The proof consists of two steps. In Sect. 2.1 we review the classification of integrable quasilinear 
systems ([3]) which may (potentially) occur as dispersionless limits of integrable soliton equations. 
In Sect. 2.2 we discuss the general procedure of the reconstruction of dispersive corrections based 
on the requirement of the inheritance of hydrodynamic reductions. This procedure is applied to 
fifth order equations in Sect. 2.3, leading to the proof of Theorem 2. 

2.1 Classification of integrable dispersionless limits 

For a system of the form 0, 

u t = tpu x + ipu y + rjWy, w x =u y , 
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the integrability conditions were obtained in [6] based on the results of [I]. They constitute an 
involutive system of second order PDEs for the coefficients (p, ip and rj, 



rj rj rj 

, ~fwWw + tyuWw - 2<PwVu + 2r)w<Pu , _ VwWu , _ Tlv)i>v 

Yuu — > tyuw — j Www 



7] 7] 7] 

u ) f]w T]u Vw 

VuU J VuW ! VwW ! 

rj 7] rj 

we assume rj ^ 0: this is equivalent to the requirement that the dispersion relation of the system 
(j3J) defines an irreducible conic. The integrability conditions are straightforward to solve. First 
of all, the equations for rj imply that, modulo translations and rescalings, one can set rj = 1, 
rj = u or 7] = e w h{u). We will consider all three possibilities case- by-case below. Notice that 
(p and are defined up to additive constants which can always be set equal to zero via the 
Galilean transformations of the initial system (J3]). Moreover, the integrability conditions are 
form-invariant under transformations of the form 

(p = <p — sty + s 2 7], ip = if} — 2sr], r) = 7], u = u, w = w + su, s = const, 

which correspond to the following transformations preserving the structure of system Q: 

x = x — sy, y = y, u = u, w = w + su. 

All our classification results are formulated modulo this equivalence. 

Case 1: rj = 1. Then the remaining equations imply ip = aw + f(u), <p = f3w + g(u), where / 
and g satisfy the linear ODEs 

r = a(f'-(3), g " = 2ag'-/3f'-f3 2 . 

The subcase a = leads to polynomial solutions of the form 

y = ju, p = {3w - ^/3(/3 + j)u 2 + 5u. (5) 

Up to equivalence transformations, the case a^O leads to exponential solutions, 

w = w + (3e u , <p = ae 2u , (6) 
where a, /3,7 are arbitrary constants. 

Case 2: rj = u. Then the remaining equations imply ip = aw + f(u), <p = (3w + g(u), where / 
and g satisfy the linear ODEs 

uf" = a(f - 0) - 2/3, ug" = 2ag' - /?/' - 1 . 

The case a ^ {0, —1,-1/2} leads to power- like solutions of the form 

if) = aw + 7 n a+1 , ip = 5u 2a+1 . (7) 
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The subcase a = leads to logarithmic solutions, 

tp = -2f3ulnu - (3u, <p = (3w + /3 2 uln 2 u + Su. (8) 
The subcase a = — 1 gives 

ip = — w + 7 In u, ip = 5/u. (9) 
Finally, the subcase a = — 1/2 gives 

ip = — -w + "yy/u, cp = Slnu. (10) 

Case 3: rj = e w h(u). Then the remaining equations imply tp = e w f(u), (p = e w g(u) where /, g 
and h satisfy the nonlinear system of ODEs, 

h" = f'-g, g"h = 2fg'-gf'-g 2 , f"h = 2hg' - 2gti + //' - fg. 

Although the structure of the general solution is this system is quite complicated, one can show 
that Case 3 cannot occur as the dispersionless limit of an integrable soliton equation. 

2.2 Reconstruction of dispersive corrections 

Given an integrable dispersionless system of the form ([3]), one has to reconstruct dispersive 
terms. This can be done by requiring that all hydrodynamic reductions of the dipersionless 
system are inherited by its dispersive counterpart [5j |6]. Following [6], we will illustrate this 
procedure with the example of the KP equation, 

u t = uu x +w y + e 2 u xxx , w x =u y . 

The dispersionless KP (dKP) equation, 

Ut = UU X + Wy, W X = Uy, 

possesses one-phase solutions of the form u = R, w = w(R) where the phase R(x,y,t) satisfies 
a pair of Hopf-type equations 

R y = fiR x , Rt = (/i 2 + R)R X ; (11) 

here [i(R) is an arbitrary function, and w' = [/,. Equivalently, one can say that Eqs. (jlip 
constitute a one-component hydrodynamic reduction of the dKP equation. Although the dKP 
equation is known to possess infinitely many A^-component reductions for arbitrary N [3 El El [TO], 
one-component reductions will be sufficient for our purposes. The main observation of [5j is that 
all one-component reductions (fTT]l can be deformed into reductions of the full KP equation by 
adding appropriate dispersive terms which are polynomial in the x-derivatives of R. Explicitly, 
one has the following formulae for the deformed one-phase solutions, 

u = R, w = w (R) + e 2 L'R XX + ^" -(fi'fjR 2 .) +0(e 4 ), (12) 



6 



notice that one can always assume that u remains undeformed modulo the Miura group [2] . The 
deformed equations (|lip take the form 



R y —fiR x 

+ e 2 {^'R xx + l -^'-^f)Rl) +0(^ 4 ), 

R t =(fi 2 + R)R X 

+ e 2 ((2/i/x' + l)R xx + (w" - ^'f + (n') 2 /2)R 2 x ) x + 0(e 4 ), 

see [5]. In other words, the KP equation can be 'decoupled' into a pair of (1 + l)-dimensional 
equations (|13p in infinitely many ways, indeed, fJ.(R) is an arbitrary function. The series in (j 12[) 
and (|13p contain even powers of e only, and do not terminate in general. 

Conversely, the requirement of the inheritance of all one-component reductions allows one 
to reconstruct dispersive terms: given the dKP equation, let us look for a third order dispersive 
extension in the form 

u t = uu x + Wy + e(...) + e 2 (...), w x = u y , (14) 

where the terms at e and e 2 are homogeneous differential polynomials in the derivatives of u and 
w of the order two and three, respectively. We require that all one-component reductions (jlip 
can be deformed accordingly, so that we have the following analogues of Eqs. (fTHj) and CE2D, 

U = R 7 w = w(R) + e(...) + e 2 (...)+0(e 3 ), (15) 

and 

R y = fj,R x + e(...) + e 2 (...) + 0(e 3 ), R t = (/i 2 + R)R X + e(...) + e 2 (...) + 0(e 3 ), (16) 



respectively. In Eqs. (|15p and ()16p . dots denote terms which are polynomial in the derivatives 
of R. Substituting Eqs. (fTB]) into (fT4"j) . and using P^j) along with the consistency conditions 
Rt y = R y t, one arrives at a complicated set of relations allowing one to uniquely reconstruct 
dispersive terms in (|14j) : not surprisingly, we obtain that all terms at e vanish, while the terms 
at e 2 result in the familiar KP equation. Moreover, one only needs to perform calculations up 
to the order e 4 to arrive at this result. It is important to emphasise that the above procedure 
is required to work for arbitrary /x: whenever one obtains a differential polynomial in [i which 
has to vanish due to the consistency conditions, all its coefficients have to be set equal to zero 
independently. Another observation is that the reconstruction procedure does not necessarily 
lead to a unique dispersive extension like in the dKP case: one and the same dispersionless 
system may possess essentially non-equivalent dispersive extensions. In particular, VN and 
modified VN equations from Theorem 1, as well as BKP and CKP equations from Theorem 2 
have coinciding dispersionless limits. 

Let us now turn to the general case of dispersionless equations of the form ([3]), 

U t = (fU x + l/jUy + 7]Wy, W X = Uy. 

The corresponding one-component reductions are of the form u = R, w = w(R) where R(x, y, t) 
satisfies a pair of Hopf-type equations 

Ry = fj,R x , R t = ((p + ipn + rjji 2 )R x ] 
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here n(R) is an arbitrary function, and w' = fi. Let us seek a third order dispersive deformation 
of system Q in the form 

u t = (pu x + ipUy + r]w y + e(...) + e 2 (...), w x = u y , 

and postulate that one-phase solutions can be deformed accordingly, 

u = R, w = w(R) + e(...) + e 2 (...) + 0(e 3 ), 

where 

R y = fxR x + e(...) + e 2 (...) + 0(e 3 ), R t = (if + # + W 2 )R X + e(...) + e 2 (...) + 0(e 3 ). 
Proceeding as outlined above we reconstruct dispersive terms. 

Remark. We point out that the formulae for dispersive deformations contain the expression 

r/w/J 3 + (ip w + f]u)^ 2 + (<A/j + 4>u)H + ¥u 

in the denominator. Since /i is assumed to be arbitrary, this expression is nonzero unless (p, if), r] 
satisfy the relations 

Vw = 0, ip w + r] u = 0, ip w + tp u = 0, <p u = 0. (17) 

These relations characterise the so-called totally linearly degenerate systems. Dispersive deforma- 
tions of such systems do not inherit hydrodynamic reductions, and require a different approach 
which is beyond the scope of this paper. 

2.3 Classification of fifth order equations 

In this Section we summarize the classification results for integrable fifth order equations 

u t = (pu x + ipUy + 7]w y + e(...) + e 2 (...) + e 3 (...) + e 4 (...), 

which are obtained by adding dispersive terms to integrable dispersionless candidates from Sect. 
2.1 (one can show that all terms at e and e 3 have to vanish). Thus, we follow the classification 
of Sect. 2.1. We concentrate on the case when the e 4 -terms contain at least one fifth order 
derivative of a or w, and skip all cases leading to third order equations which were already 
classified in [B]. 

Case 1: We verified that the exponential solutions Q do not survive, so that all non-trivial 
examples come from the polynomial case (|5|), 

77 = 1, ip = ju, f = /3w - ^/3(/3 + j)u 2 + 5u. 

A detailed analysis of dispersive deformations leads to the constraints 7 = j3, 5 = 0. Modulo 
rescalings, this gives BKP/CKP equations. 

Case 2: One can prove that none of the logarithmic cases ([8]), ([9|) and (|10p survive, so that all 
non-trivial examples come from the power case ([7]), 

77 = u, ip = aw + ju a+1 , (p = 5u 2a+1 . 

The further analysis leads to the only possibility a = —3. Modulo rescalings, the case 5 = 7 = 
gives the HD§ equation. The case of nonzero 5 and 7 leads to the new equations E3 and E±. 

Case 3: One can show that no examples from this class possess fifth order dispersive extensions. 
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3 Lax pairs 

For the reader's convenience, in this section we bring together Lax pairs for all equations ap- 
pearing in Theorems 1-2. We emphasise that our classification scheme does not assume the 
existence of a Lax pair: these come as the result of direct calculations once the classification 
is completed. We refer to [T7] for an alternative approach to the classification of integrable 
systems in 2 + 1 dimensions based on postulating the structure of a Lax pair. 



3.1 Third order equations 

Since both KP and modified KP equations are particular cases of the Gardner equation, we will 
skip the first two examples. 
The Gardner equation, 



/3 2 

u t = (fiw - —u l + Su)u x + w y +e 2 u xxx . 



.2 



possesses the Lax pair (TTJ 



6 S 

e 2 ip xx + —=(if) y - f3uip x ) + -uip = 0, 
V3 6 



etpt = ^ipyyy - V3(3e 2 (2-p xx + u x ip x ) + e((3w + —u 2 + Su)ip x + e-u x - J^ - " 2 + ~^J^ W - 
The VN equation, 



U t = (uw)y+e 2 Uyyy, 



possesses the Lax pair [161 HE] 

1 

3 



<?"Pxy + 7;U1p = 0, 



The modified VN equation, 



3u 2 



U t = {uw)y+e [Uyy---^] , 



possesses the Lax pair [T] 

2u Px + 3 



4 u 



2 2 V 

e Wxy - e t^^x + ^uip = 0, 



1p t = e Ipyyy + Wt/jy + -Wyljj . 

The HD equation, 

Ut = —2WUy + UWy ( — 

possesses the Lax pair (TTJ 



4e 2 , / 2V3ew 6e 2 u x \ , 

ipt = —»Vxxx + 5 1— Wxx- 

u° \ u A u I 



9 



The deformed HD equation, 



5 e 2 (\ 

Ut = —U X — ZWUy + UWy [ — 



w 

possesses the Lax pair [6] 

£ 2 e 5 

4e 2 , ( 2V3ew 6e 2 u x \ , S , ( 35u x y/35w" 
* = + ^ + ^ + 1^-1^- + 1^2" ^ 

The Equation Ei, 

Ut = (/3w + ftu 2 )u x - 3j3uu y + w y + e 2 [B 3 (u) - /3u x B 2 (u)], 
possesses the Lax pair [6] 

e 2 V^ = e 2 /3wp xx + i^, 

■0 t = e 2 ftu 3 ip xxx - e 2 ip yyy + 3e 2 ''ft' l uu y ip > xx + j3wip x . 

The Equation E 2 , 

n t = ^ftu 3 u x + (w - 3/3u 2 )u y + uw^ + e 2 [5 3 (u) - f3u x B 2 {u)], 
possesses the Lax pair [B] 

1 



e 2 ipxy = e 2 (3uil) xx + ^utp, 



, , ft 

1p t = e /9 U 1p XXX - € Ipyyy + 3C ft UUy^ XX + — U l\) X + Ufl/fy + fiUUylp . 



.2o3„,3„/, 2„/, i nJo2„,„, „/, , K „,3 



3.2 Fifth order equations 
The BKP equation, 



u t = 5(u 2 + w)u x + 5uu y - 5^ + 5e 2 (uu xxx + u xxy + u x u xx ) + e A u x 



possesses the Lax pair [TT] 



tfiy + utp x + e = 0, 

tpt + 5(u 2 - w)t/j x + e 2 (15ut/j xxx + 15u x ip xx + lOtt^s) + 9e 4 VWxx = 0. 



The CKP equation, 



5 

u t = 5(u 2 + w)u x + buuy - 5w y + 5e 2 (uu xxx + u xxy + -u x u xx ) + e 4 u x 
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possesses the Lax pair 



ip y + uip x + -u x ip + e ip xxx = 0, 

1 



45 



35 



tp t + 5 (it - w)ip x + 5(uu x - -u y )ip + e (15m/w + —u x ^ xx + —u xx ip x + 5u xxx i^) + 9e tp. 



The HD equation, 



ut = 15wu v — 5uw v + 5e z 



^xxy 3 / U x Zly 



U 2 U V U 2 J x 



e 

2U 2 



possesses the Lax pair [IT] 

e 2 



2 1" 4 9 15/1 

% + 15e nlpxXX 4" £ zIpXXXXX ~ 45 aYXXXX H o 1 ~~ 7 



. , l^xxx 



= 0. 



The Equation E 3 , 

Ut =4 r Y 2 —^ + 5(3iu — \)u v — 5uw v 



4-56 2 

possesses the Lax pair 



J_ / J_\ + Itzzy _ 3 /% 



2u 2 \u 2 



xxx 



l x Uy 



u \ U A 



2u 2 In 2 , 



7 e 
^x + -^ripxxx = 



(67 2 15710 "\ / iu 7 \ 3^yit ^7 Z' ^ 

—5 FT" ^ + 15e 2 ( r ) 1p xxx H ^Vxx + -fT ( — ) ^rr 
u° u a J [Vr uv u° V It / x 



+e 4 

The Equation E 



AK u x 15 / 1 \ 

"WxXXXX 45 "F^XXXX ~^ o I 9 I WxXX 
U U \ U J XX 







K 




15 




u 3 


U 2 



itt =47 2 — | + 5 (3iu — -tt W — 5uw v 
u° u z 

7 ( 1 



-5e 2 



7 - 7 1 u 



3u \« 3 /xxx n? 



^ "U X Uy " 



— e 



■ xxy 

1 fl\ _ 15 / (2uu^-3n 2 ) 2 
16 V u 8 



U A J x 



possesses the Lax 



2u 2 
pair 



/7it x u y \ 7 



u* 5 



3 tt x:r 9 u 2 \ / 1 / u xx \ 3 -u 3 \ 



11 



+ 



15(u 3 u y + 2ju x )w lO^Uy lQ^ 2 u x 5 



2u A 
lh(wu 2 — 7) 



+ 



( 607^ 15u 



6 



U7 y V + I 6 



?7," 



15^ ] V, + 



,6 



+ 



2v? 



+ 
+e 
+ 
+ 



/ 307 /M^X 

57 {-WJ 



I5u x u y (90uu x 



135<)w 



\^x 



2u A 

157 ^ V? 

x U 

45u 



+ 



u 

45 fu xxx \ { u x ii xx \ 

/ 1 \ 195 /u x u xxx \ 135 /i&A 4605 /u 2 

~ W„ + "2" \—a^)x + —\lf) x - — \- 



u 



6 



n u J . 2 

tpxxxx + 



4u 5 



u>xy 
U 3 / x 



1px + 



15 

+ 2T 



1 k)Uxx 

105n| 
2u 7 



+ 



225< 



45ni! lbu'iuy , 
+ -r^r^ H ^4^- 1 V> 



4m 6 
r ( 30 



u xx 
u 6 



2u 5 



+ 180+ | 



+ 



1pxx + 



825 /«■ 



-I -± 

6w 4 V 



235u 4 ' , 



+ 



3375 /m 4 



4 x 3645u^ 
+ 



165^4^ 
2u 8 



j; 



8w 10 



We do not exclude a possibility that simpler Lax pair can be found in this case. 
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